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NOMENCLATURE 
Average sidewall wetted area, 
starboard side 
Average Sidewall wetted area, port side 


Added mass coefficient in roll force 
equation 


Added mass coefficient in pitch force 
equation 


Average wetted sidewall area of the bow 
Average wetted sidewall area of the stern 


Added mass coefficient in yaw force 
equation 


Sideslip angle 
Coefficient drag on x-direction 
Coefficient drag in y-direction 


Sidewall roll moment lumped parameter 
coefficient 


Bow pitch force lumped parameter 
coefficient 


Sidewall roll force 

Sidewall starboard buoyancy force 
Sidewall port buoyancy force 
Stern buoyancy force 

Bow seal pitching force 

Bow buoyancy force 

Gravitational acceleration 

Moment of inertia about x-axis 


Moment of inertia about y-axis 
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fete 


ie ele; 
ft-slug 


ft 


fea 


slug 5° 


li 
ft = 


rad 


2 2 
1b;s lb/ft 
2 2 
lb;s lb/ft 


lb,s*/ft* 
non- 
dimensional 

lbs 

lbs 

lbs 

lbs 

lbs 

lbs 

Ft/s* 
slug us 


slug ped 
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Moment of inertia about z-axis 
Track deviation performance index 
Control effort performance index 


Overall performance index 


Summation of moments about x-axis 
Length of sidewall 
Actual draft of port sidewall 


x-direction displacement of hull 
drag centroid 


Average draft of SES sidewall 

Average draft of bow seal 

Length from center of gravity to stern 
Average draft of sidewall 


Pitch moment lever arm for bow sidewall 
buoyancy force 


Pitch moment lever arm for bow seal force 
Mass of the rigid ship 

Summation of moments about y-axis 
Summation of moments about z-axis 

Plenum pressure 

Lumped drag centroid point 


Lumped drag centroid point 


Roll acceleration 
Pitch rate 

Pitch acceleration 
Yaw rate 


Yaw acceleration 


11 


Slug gee 


Eeomeee 


Pace sec 


non- 
dimensional 


lbs-ft 
Ee 
Et 


ete 


fit 
Bit 


Pit 


ae 


ft 

fae 
lbs-ft 
lbs-ft 
ibs ft" 
non- 
dimensional 
non- 
dimensional 
pay eee 
rad/sec 
Pa eee. 


rad/sec 


pad sea- 





Fforw 
eae 


yaw 


Turning moment lever arm of no. 1 engine 
Turning moment lever arm of no. 2 engine 
Turning moment lever arm of no. 3 engine 
Turning moment lever arm of no. 4 engine 
Total thrust magnitude on no. 1 engine 
Total thrust magnitude on no. 2 engine 
Total thrust magnitude on no. 3 engine 
Total thrust mangitude on no. 4 engine 
Total forward thrust vector of effectors 
Total side thrust vector of effectors 


Total turning moment generated by 
effectors 


Velocity in x-direction (surge) 
Acceleration in x-direction 
Velocity in y-direction (sway) 
Acceleration in y-direction 

Total velocity 

Width of bow seal 

Summation of forces in x-direction 
X coordinate of SES 


NAV 


XA Velocity Of SES 


Summation of forces in y-direction 
Ywav coordinates of SES 


YNWAV velocity of SES 
Effector angle commanded 
PoaaSetor angle Of no. lL nozzle 


Effector angle of no. 2 nozzle 


1 es 


ft 
fete 
i 
Et 
lbs 
lbs 
lbs 
lbs 
lbs 
lbs 


lbs ft 


ft/sec 
ft/sec 
ft/sec 
ft/sec 
ft/sec 
ee 
lbs 
ft 
Ete 
lbs 
ft 
ete 
eel 
rad 


rad 
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EEteceon angle Gf no. 3 nozzle 
Effector angle of no. 4 nozzle 
Heading angle of SES 

Roll angle of SES 

Pitch angle of SES 


Density of sea water 
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rad 
rad 
rad 
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I. INTRODUCTION 


Extensive analysis of the dynamics of the Captured Air 
Bubble (CAB) Surface Effect Ship (SES) has been accomplished 
in the past. Such analysis has produced equations of motion 
of the SES type vehicle in six degrees of freedom [Ref. 1]. 
Analysis at the Naval Postgraduate School has centered 
around the XR-3 and 100-B SES [Refs. 2,3], with current 
analytical effort being directed toward the 3K-SES [Ref. 4]. 
The primary tool for study of the 3K-SES characteristics 
is the highly comprehensive five degree of freedom (DOF) 

SES DBS Program developed by the David W. Taylor Naval Ship. 
Research and Development Center [Ref. 5]. <A version of 
this program used at NPS is called DBSIM5D. This data base 
program allows very detailed study to be conducted on the 
3K-SES response to precomputed control inputs such as step 
or ramp thruster commands and their associated accidental 
failure possibilities. 

A natural extension in the study of control of the 
3K-SES is to evaluate the vehicle's response to an on-going 
"man generated" control input on a real time basis. The 
rationale for developing a "man-in-the-loop" SES simulator 
would be to 1) test the man in responding to accidental 
failure situations that may be experienced during opera- 
tional conditions, 2) train personnel in operator control 


techniques unique to SES type vehicles, 3) provide a means 
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to test hardware to be utilized in SES craft, and 4) provide 
a design tool for development of control schemes or system 
structural changes to improve performance. 

Accordingly a simplified, non-linear five degree of 
freedom Real Time Simulator (RTS5D) for the 3K-SES was 
developed such that a continuous, on-going, observable, 
real-time solution is interfaced with man-generated control 


inputs. 
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LI. EQUATIONS OF MOTION 


A. COORDINATE SYSTEMS AND ASSUMPTIONS 
It was recognized that the RTS5D would require accurate 
and yet simplified equations of motion when compared to 
those of the more extensive DBSIM5D. The rationale for 
requiring simplification was to start with a right hand 
Side (RHS) that could be utilized effectively in a real 
time simulation. To satisfy this requirement a model was 
selected that would use point source of drag forces con- 
Sisting of only a few lumped-coefficient terms. The 
equations used are based upon those developed in Ref. l 
which follow this concept. Additional assumptions used in 
the simplified 5 DOF equations were as follows: 
1) All accelerations are measured at the center of 
GEavi ty. 
2) Vehicle is "free-to-heave". 
3) All cross coupled moments of inertia are zero. 
4) All cross products of angular velocities in force 
equations are zero. 
5) All roll and pitch angles used in force equations 
are subject to small angle approximations. 
6) Calm water conditions. 
7) Mass and mass distribution of vehicle are constant. 


8) Effect of changes of aerodynamic force are neglected. 
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Utilizing the above assumptions, the simplified 5 DOF 
equations of motion are developed using the coordinate 
systems shown in Fig. 1 and Fig. 2, the force vector 
diagrams shown in Fig. 3, Fig. 4A, Fig. 4B, and the 


Following equations from Newton's laws of motion: 


SURGE m(u - vr) = xX 

SWAY Tey meta), 

YAW eae = oN 

PITCH I, q = M 

ROLL I p nr 

where: 
m = mass of the rigid ship slugs 
u = velocity in x-direction (SURGE) pee See 
Vv = velocity in y-direction (SWAY) Rt e/ sec 
i = angular velocity about the z-axis , deg/sec 
p = angular velocity about the x-axis , deg/sec 
ay = angular velocity about the y-axis , aeg/sec 
Z 7 moment of inertia about the z-axis , lb -£t* 

Iy = moment of inertia about the y-axis, lb-ft" 
Ty = moment of inertia about the x-axis , lb_-£t* 
X = summation of forces in x-direction , 1b, 
y = summation of forces in y-direction , 1b, 
N = summation of moments about the 


Z-axis ‘ ft-lb. 
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M = summation of moments about the 
y-axis Ae ee Lb 


K = summation of moments about the 
X-axls , f£t-lb 


Additionally the following navigation relationships 


are utilized (see Fig. 1). 


X = u cos J - v sin yp 
XY, = 4 Sin »y + v cos vp 
VV = 7 + ye 
Ss 
-l 
8 = tan (-v/u) 
where: Ww = heading angle 
B = drift angle 
6 = thrust vector angle 
Also, i f eit een cit 
oS f ¥, dt if arch 
meee | rode i eed t 
Eaeea) pdt ¢ = fpdt 
q = f{qdt e = fq dt 
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NAV 


NOTE: 1) w=0 when vehicle heads parallel to + AN AY 
axa S. 


2) Vehicle in figure above shown with negative 
thrust vector angle 6, positive w, negative 
B, negative sway velocity v, 1.e., ina 
eens eU Tl. 


Figure l 
Definition Of Coordinate System (Part I) 
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tz, 
EG, / 
as 
u 
7 
Tae 
/ 
/ +p 
/ 
4 | JL, *q 
sy fy 
Z, 
-V +o ey 
a ae ea 
—ne; 
NOTE: Direction for positive velocity and 


angular rates are shown 


Figure 2 
Definition Of Coordinate System (Park IT) 
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BOR L 
SIDEWALL 


STARBOARD 
oS IDEWALL 


a PN | > J 
7 a ae acl. lo ala 
T, C085) 4 


NOTE:1) p and p' are equivalent point force centroids. 
ay. 6 6 


79 99> Sg» 85 9 are thrust vector angles and 
are not a a to be equal. 
Sipe s.. L are thrust magnitudes and are 


Oe hatha +8 Percale le. 
4) Force directions shown are for a right turn 


Figure 3 
Surface Effect Ship (Top View) 
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Sp sw 


Pert iMelUrieronces SACtIng, 1.6€., 64,5 6 


7 Se 


S49 are negative 


TorpE = ee ee oe ign) ta 


Figure 4A 
Surtace Effect Ship (f&tern View) 


Zz 





+Z 





MOreE: wal) Tropwrt7 cos 57 tT gcos SgtTgcos 54tT 


(see Fig. 3) 


cos 6, 


1 0 


Figure 4B 
SimpaceseGiiectsShap (Side View) 
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B. FORCES AND MOMENTS 
l. Surge Forces (See Fig. 4B) 
The forward acceleration was determined to be 
generated by the summation of thrust vectored in. the forward 


direction T The resultant acceleration of the SES was 


forw- 
assumed to be opposed by a retarding force exhibiting a 
"“velocity-~squared" characteristic. Additionally, Newton's 
laws of motion specified a "v-r" product that contributed 

to the drag component of the equation when the vehicle was 
in a turn (centrifugal force reaction term). The simplified 
equation of motion for the surge acceleration follows from 

a summation of forces in the u direction as depicted in 

Fig. 4B. 


ue ( 7) ee (C, a” /m) 1 Ns 


T forw 
2. Sway Forces (See Fig. 4A) 

The side acceleration of the vehicle was analyzed 
to be the result of the summation of the thrust vectored 
parallel to the vehicle's y-axis and a retardation "velocity 
Squared" phenomena such as that produced by the cross flow 
drag term of the sidewall as described in Ref. 6. This 
retardation force is augmented by a "u-r" product as 


specified by Newton's Law. 


fit vivi|/m - ur 


Ge 
hy 


Side 
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Note the requirement to force the Cv /m term to 
maintain the sign of the sway velocity such that both left 
and right turns may be computed. 

3. Yaw Moments (See Fig. 3) 

The yaw acceleration was a turning moment summation 
which was defined to be influenced by three primary forces 
operating over moment arms. It is of interest to note that 
the third term is an "added mass" term which is considered 


to operate on the same pressure point as the C viv tee iain. 


DY 
The term va is a summation of moments generated by the 


four thrusters as shown in Fig. 3. 


a. = s,T. cos 6 tS. 1 GOS) on — S-1. Cos 6 


yaw 4° 7 3°8 8 1°79 9 
_ Soli g cos 836 = & oT. Sin 4 - Lot, Sin 52 
- & oT, sin 5 = & TI 6 Sin Si9 

eee SW Taw tz + A,,uvd, /1, 


4. Pitch Moments (See Fig. 4B) 
The pitch acceleration was assumed to be the summa- 


tion of moments generated by forward thrust, T the 


forw’ 
buoyancy of the sidewalls of the SES, Fy and Fos and a 
vertical force generated at the bow of the vehicle, Fo: 


This vertical bow force was defined to be a lumped parameter 
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term which modeled the reaction force due to plenum pressure 


acting against the bow seal. 


BL 


average draft of bow seal 


Ld a (C2, tan 8) /2) 2, 


ale 
LO, - ((4, tan 6)/2)2, 
Ay “x, 0g 
Aw LX. 0g 
CapP We (%q - 23, tan 6) 


average draft of sidewall 


plenum pressure 


average wetted sidewall area of the stern 
width of one sidewall 

width of bow seal 

average wetted Sidewall area of the bow 
added mass coefficient 

water density 

B25 ee 


lever arm of bow seal 
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2 
il 


3 lever arm of buoyancy force 


Qe 
| 
4 
rh 
O 
KX 
= 
<s 
oS 
UJ 
O) 
bt 
NO 
WwW 
Oo 
rs 
N 


- Fi 2. - A,3uq)/1y 


Note the added mass term A,,uq which was required 
to provide damping to the pitch moment and is specified 
iaeRet. lt. 

It was found that the response of the SES toa 
pitch perturbation without the added mass component was 
undamped and approximately sinusoidal. Unlike the yaw 
equation where the added mass term was a small contributor 
to damping, the pitch added mass term was found to be of 
Significant importance in modeling the known pitch motion. 

5. Roll Moments (See Fig. 4A) 

The roll acceleration equation was analyzed to be a 
summation of moments generated by buoyancy forces of the 
port and starboard sidewalls, ESD and Bae! the thrust 


vectored parallel to the y-axis ofthe vehicle, T a Side 


Sage. 


force C,, v|v|, and a lumped coefficient vertical force, 


: The vertical force eae was defined as the force 
generated in a turn due to the sidewall curvature (dead 


rise angle) acting against the cross flow of water. 
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Psp ~ 9 Ayp*dp 
ss ~ 9 Ars * as 
Few = Cpgp VII 
“ay = average draft of SES sidewall 
ap = “ay = ee tan 
Las = “a, + Brees tan 6 
Arp = average wetted area port 
as = average wetted area starboard 
p ~ (UF oo - P ss) sw - side ae * “py we *y 


s P sw * Sw 7 ai aoa 
Again note the added mass term A,,uUp/I, which was 
found to be essential in modeling the damping phenomena of 
the roll motion. 
In summary, the simplified 5 DOF equations of motion 


used in RTS5D for the 3K TON SES are: 


/m=-C u* /m + vr [Fig. 4B] 


SURGE u DX 


T 


forw 


SWAY We Tosge/™ ee vivl/m - ur [Fig. 4A] 
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YAW = Taw’ =z + Coy viv 2/1, + Azouvd, /I, [Fig. 3] 
E 2 
= = y 
PITCH q Gi a ee eee Uz 
~ Fi 3 = A,,uq)/I, LPig=. 243) 
ROLL p = nee 6 ome eee i aac, viv, 


~ Be as — A, up)/I, [Fig. 4A] 


D. PARAMETER IDENTIFICATION 

These equations are an extension of the 3 DOF flat turn 
SES model developed by Gerba and Thaler in Ref. 4. The 
DX” ~ DY’ end W is 


described in Ref. 4 and repeated here for completeness. 


identification of craft parameters C C 


The surge drag coefficient C is determined by selecting 


DX 


a steady state turn condition, where T lie Va8 oe and 1 


forw’ 


are known. 





2 
- r Cue) 
= LOEW DX nen. 
m m 
from which 
TE + mvr 
- = forw 
DX 2 
u 
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The sway drag coefficient Coy is determined by the same 


steady state turn condition which requires 


ss co. Cc 
Vv = 0 _side _ bY viv = ur 
m m 
where 
C ae eae 
DY 
viv | 


The sway drag moment arm follows utilizing the yaw 


acceleration equation 


mo. . ees . Coy VIVI2, . A, uv 2 
T a wares St. 
Z Ty Ty 


DY vij[v| + A,juv 


a2 


The addition of the pitch and roll equations introduces 


DP ZP° These are 


eaSily solved using known constants ae ae Coy? Loe een 


additional lumped coefficients C and Cy 


Wa! Pp Xs and steady state values of Psp’ Fog! Fie Pos 


Au ren L 


lpeveednd 2s. [Ft 1S Significant to note that 
forw 


Side’ 
in a steady state condition the angular roll rate, p, and 


angular pitch rate, q, are both required to equal zero. 
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This is in contrast to the yaw rate, r, where a finite 
steady state value is desired in a turn. Thus the added 


mass terms of A, 3uq and A,,up in the pitch and roll equations 


ou 
are not utilized in the determination of craft parameters; 
their function is strictly confined to damping of their 


respective accelerations. Therefore, utilizing the pitch 


equation with known steady state turn values, 


coy pt “ppepte <a “31 “am 8) *3] 
+ F.2, - C op = P-¢.)/71 
Z2°3 DX 2 es Y 
yields 
Cus F.g. - F.t. - T Q 
‘ Beeox 2)" 173 E38 forw “p 
MMi (0 1- C.. tan oj... 
DP PL. Ls Lot tan O)2%34 


The roll equation under steady state turn conditions 


yields 

= = ~ = . + 

P : (Pep Pos) tw "side t “py wy ty 

- pgp Vivi 2,,,)/T, 
from which 
= ~ - : X 
Cpzp (Foo Fos) iw T side ae u py mails | 2)/ 
Wee es) dow? 


ah 





where 


FE 
SP 


ss 


side 


port sidewall buoyancy in pounds [Fig. 4A} 
starboard sidewall buoyancy in pounds [Fig. 4A] 


thrust vectored parallel to vehicles y-axis 
i pounds (Fig. 3] 


Sway steady state velocity in feet/sec 

sway drag coefficient for steady state condition 
Sway drag moment arm in feet PEGs 

effector thrust moment arm in feet [Fig. 4A] 
Surge drag moment arm in feet (Eagee. 4A] 


sidewall buoyancy force moment arm in feet 
[Fig. 4A] 
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III. INTEGRATION METHOD 


A. ERROR ANALYSIS 

A rectangular (RECT) type integration was used to 
compute the velocity and positional variables of the force 
equations. Selection of the type of integration over a 
more complex type such as Runga Kutta, Variable Step enabled 
the loop delay time due to integration to be held to a 
minimum. 

An analysis of the error introduced into the solution 
of the simplified 5 DOF equations through the use of 
rectangular integration vice Runga-Kutta, Variable Step was 
made. The simplified 5 DOF equations of motion were 
programmed on the NPS IBM 360/67 computer using the 
Continuous System Modeling Program (CSMP). One series of 
computer runs was made utilizing rectangular integration 
over a Spectrum of step sizes ranging from 0.301 secs to 
0.25 secs. A second computer run was made utilizing Runga- 
Kutta, Variable Step type integration. The steady state 
results are shown in Table I. Examination of the final 
value magnitudes established that no Significant steady 
state error was due to rectangular integration for the time 
increments utilized. 

Transient error of the variables, r, $6, and 9 was 


examined. A series of computer runs uSing incrementally 
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larger rectangular integration step sizes was compared 
against an identical run uSing Runga-Kutta Variable Step 
type integration. The first peak absolute magnitude 
results are shown in Table II. The RTS5D had an iteration 
time of .084 seconds therefore the "transient" magnitudes 
were in error by approximately 2% to 3% due to rectangular 


integration. 
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TABLE II. 


TRANSIENT FIRST PEAK ERROR ANALYSIS 
FOR RECTANGULAR INTEGRATION 


DELT 16 b Q 
.050 1.42% 1.64% 2.58% 
Ov SD Vs 6 a8 2.413% 2.84% 
2EOO Zi Los Boos 3425% 
er25 2.353% 3.65% 8.57% 
Note: 1) Step effector angle input of 25° 


and velocity (u) of 56 knts. 


2) All percentages are referenced to 
Runga-Kutta variable step magnitudes. 


36 





IV. IMPLEMENTATION 


A. PHILOSOPHY 

Given the decision to develop a real time motion analysis 
Simulator for the 3K-SES, the various options available 
were examined. Prime requirements for such a simulator 
were the necessity for real clock interfacing, suitable 
input-output ports, and a dedicated Central Processing 
Goat (CPU). 

The “all analog" computer simulation was particularly 
attractive from the aspect of interfacing man generated 
control inputs into an on-going real time solution. With 
the analog approach instantaneous solutions are obtained 
from patchboard amplifiers. Two considerations which did not 
lend themselves well to the all analog technique were how to 
Maintain high accuracy solutions in wide ranging variables 
when constrained to a voltage limit of plus or minus 100 
volts, and how to output the system variables in such a 
manner that a man could physically interpret what was 
happening and respond with the required corrective, timely 
control effort. Additional supplementary problems associ- 
ated with the "quarter-square" multiplication error and 
amplifier drift/run duration limitations gave rise to 
examining other possible approaches. 

The "all digital" computer simulation offered high 


accuracy with wide ranging variables and a capability to 
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interface with two graphics display consoles. Additionally 
a line printer was available to provide "hard copy" of 
variables of interest as the solution was being processed. 
Two major disadvantages were the serial and discrete nature 
of the integration method of a digital computer, and the 
lack of high voltage, high accuracy potentiometer controls 
to allow the man to influence the motion of the vehicle as 
the solution was being presented to him. 

A desire to take advantage of the best features of both 
the analog and digital alternatives led to a decision to 
develop the real time simulator (RTS5D) of the 3K-SES 
utilizing the simplified 5 DOF equations of motion developed 


in Section II. 


B. REQUIREMENTS 
The following criteria were required of the RTS5SD: 
1) The simplified equations of motion would be solved 
and the results output on areal time basis. 
2) The solution would be subject to real time control 
efforts generated by a man observing the output. 
Development of the simplified 5 DOF equations of motion 
and the use of rectangular type integration for simulating 
the 3K-SES ensured that loop delay attributable to digital 
computational needs would be minimal and thus the solution 
would be compatible to a real time environment. 
The requirement to output the results on a real time 


basis such that a man observed and responded with real time 
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control efforts necessitated the use of graphics terminals 
interfaced with a digital computer. Large volumes of data 
and graphical cues were presented to the pilot and updated 
on a real time basis in much the same manner as found on 

the bridge consoles of the 3K-SES shown in Ref. 7. Software 
algorithms designed to present "alpha-numeric" data and 
graphic perspectives corresponding to the vehicles current 
State variables enabled the observer to respond on a real 
time basis. 

Real time man generated control requirements necessitated 
the need for a hybrid environment and its associated inter- 
face between a digital computer and an analog computer. 

The RTSSD utilizes four (4) computer systems: One (l) 
medium size general purpose digital computer (XDS-9300), 
two (2) small graphics digital computers (AGT-10), and one 
(1) analog computer (CI5000). The XDS-9300 has a CPU 
capable of handling 32K words of main (core) memory with a 
magnetic drum for secondary storage. The XDS-9300 may be 
programmed in assembly language or FORTRAN. The CI-5000 

is an analog computer with a full complement of digital 
logic. Hardwired trunk lines linking the XDS-9300 and 
CI-5000 allow hybrid computer operations. The two AGT-10's 
are small general purpose digital computers with 8K of main 
memory and magnetic disk for secondary storage. The AGT-10 
systems have the capability of accepting input and of 


displaying output on their CRT screens. Both the CI-5000 
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and the AGT-10 have the capability of communicating with 
the XDS-9300 directly, but cannot communicate with each 


other except via the XDS-9300 [Ref. 8]. 


C. HARDWARE DESCRIPTION 

Hardware layout for the RTS5SD is shown in Figure 5A. 
A "joystick" type high voltage potentiometer attached to 
the pilot's seat was hardwired to the CI-5000. This joy- 
stick served the purpose of the pilot's effector angle 
controller. Left or right stick transmitted an appropriate 
nozzle deflection angle in the form of a voltage level to 
the CI-5000 where it was quantized in an analog to digital 
process and introduced into the XDS~-9300's ongoing solution 
of the simplified 5 DOF equations. In like manner the 
voltage levels of six other "dial" potentiometers used to 
control thruster output levels and program modes were 
introduced into the solution process as deemed appropriate 
by the operator. The six potentiometers were Weeatad en a 
"thruster console" (See Fig. 5B) within reach of the pilot 
(See Fig. 5A). Additional hardware included a standard 
digital computer line printer, an analog 8 channel strip 


chart recorder, and an analog X-Y plotter. 


D. SOFTWARE DESCRIPTION 
if Introduction 
The RTSSD software package was a FORTRAN IV digital 


computer program which consisted of a main program and two 
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major subroutines. The program gave the operator the option 
to start, stop, and hold the dynamics solution process. 

The software package controlled a multiplexing algorithm 
used on the graphics display consoles and sampled all 

analog control and timing inputs. A program flow chart is 
shown in Fig. 6A and a detailed flow chart for the graphics 
multiplexing is shown in Fig. 6B. 

The program began by initializing all program 
variables and parameters and then requested from the operator 
any parameter changes. The request and input of changes 
were made on a TTY-35 Teletypewriter (See Fig. 5A). The 
program then started the real time 1000 Hz clock timing 
Circuits in the CI-5000 analog computer and entered the main 
program iteration loop. In the main loop the sampling of 
analog control inputs by the pilot was accomplished by 
subroutine ADDA. Craft force noua ous were computed and 
integrated. A perspective view of the simplified form of 
the SES and an electronic road were utilized in the RTS5D 
to give real time visual cues to the pilot on the status 
of the vehicle's states. This section of the display gave 
the pilot his primary guidance information and was updated 
on the pilot's graphics display console on every iteration 
of the main program loop. The program then evaluated 
whether an accidental failure mode had been selected by the 
Operator. If the failure mode was selected the program 


searched the craft's states to see if the recovery criteria 
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had been satisfied through the maneuvering inputs of the 
pilot. If the recovery criteria was fulfilled, the program 
went into a hold status and awaited further instructions 

by the operator on when to restart the run. If the recovery 
criteria had not been fulfilled the program continued the 
main loop by incrementing a graphics multiplex index and 
entering the graphics display multiplex algorithm to search 
for and display the appropriate information array associated 
with the current value of the graphics index I (See Fig. 6B). 
A graphics multiplexing algorithm was established to reduce 
the iteration time required for the main loop. Each informa- 
tion array contained approximately 50 elements associated 
with positions on the display screens. Each array had a 
display time which was required to encode and project the 
array on the screen. This display time varied directly 

with the number of words in the array. The prioritization 
of the graphic display console information arrays was made 
such that the perspective of the SES and the electronic 

road were displayed on each iteration and were not subject 
to the program's graphics display multiplexing routine. 

All other information arrays were multiplexed on a time 
Share basis such that on any main loop only two information 
arrays were projected on the screen, the SES/road and the 
array designated by the graphics display index I. This 
concept was mandatory due to the large volume of data/ 


graphics being presented. It is instructional to recall 
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the fact that all variables were computed and updated in 
Ehe memory of the XDS-9300 on every loop, whereas the 
display routine was multiplexed on a time share basis. 
The main loop iteration time was 84 msec. Information 
arrays in the multiplex routine were updated on the graphics 
display consoles every 2.0 sec. On exit from the graphics 
display multiplex algorithm the program read the value of 
the real time clock and then reset the clock to zero. From 
this the A time magnitude was determined. Following this 
the program searches for print, restart, and hold flags 
which are set by the operator. The program loops back to 
the AD/DA subroutine and begins a new iteration. 

2. Graphics Techniques 

The two AGT-10 Graphics Display Consoles were 
utilized to display real time data to the operator such that 
he could control the craft. The process of transmitting and 
displaying desired data to an AGT-10 was accomplished 
utilizing FORTRAN callable graphics subroutines in the 
XDS-9300. One graphics subroutine TEXTO was used to transmit 
alpha-numeric data and one graphics subroutine GRAPHO was 
used to transmit graphical data. 

Subroutine TEXTO transmitted alpha-numeric data to 
the AGT-10 screen. The positioning of the desired data on 
the screen itself was accomplished by utilizing a TEXTO 
screen coordinate system of dividing the 12 inch by 12 inch 
AGT-10 screen into 40 equally spaced lines from top to 


bottom with each line having room for 96 individual 
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characters on it as shown in Fig. 7A. Subroutine TEXTO 
statements analogous to FORTRAN "WRITE" and "FORMAT" state- 
ments were used to specify the letter character or variable 
Magnitude to be displayed and their position on the display 
screen. 

Subroutine GRAPHO was used to transmit the "picture" 
portion of the output to the display screen. The coordinate 
system utilized by GRAPHO is a normalized X-Y system with 
the ordinate and abscissa intersecting in the center of the 
screen as shown in Fig. 7B. A point on the GRAPHO coordinate 
system corresponding to (1.,0.) would be located at the right 
edge of the screen centered between the top and bottom of 
the screen. A line could be generated on the screen by 
specifying the line start coordinates, line stop coordinates 
and commanding a "draw". The AGT-10 Weald connect the two 
designated points and a line would be displayed on the screen. 
Concurrently if only two point sources of light had been 
desired to be displayed, the command "jump" could be 
utilized such that the AGT-10 would neeeeomece the second 
coordinate point with the first coordinate point. In this 
Manner a series of commands was generated via subroutine 
GRAPHO whereby a picture was drawn on the AGT-10 screen 
using an array of sequential coordinate points and their 
respective "draw" or "jump" instruction. 

3. Analyst Graphic Display Console 
The display presented to the analyst during the 


conduct of a run is shown in Fig. 8. The magnitudes of the 
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Figure 7B 
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the variables were computed by the XDS-9300 on every iteration 
of the program. The XDS-9300 then followed a display multi- 
plexing algorithm (See Fig. 6B) which updated the variables 
every 2. seconds. The variables were categorized by position, 
velocity, and acceleration. Each category had three lines 
of data displayed. The "present" data line displayed the 
real time magnitude of the variable column it was located 
in. The “maximum" and "minimum" were the max and min 
magnitudes that had been computed for the variable during 
that particular run. At the start of a new run all variables 
are initialized (See Fig. 6A). 
4. Pilot Graphics Display Console 

The display presented to the pilot during the conduct 
of a run is shown in Fig. 9A and Fig. 9B. The pilot display 
contained three major sections, VARIABLE DYNAMICS, SEA TRACK/ 
EMERGENCY PERFORMANCE INDEX, CRAFT PERSPECTIVE/NORTH ROAD/ 
NAV DATA. 

a. Variable Dynamics 

The time histories of the accelerations and their 

respective velocities were presented in graphical form. An 
acceleration plot and its respective velocity plot shared 
the same electronic strip chart such that one was super- 
imposed over the other. The acceleration plots were 
displayed as "light" toned lines whereas the velocity plots 
were displayed as "heavy" toned lines. As the RTS5D program 


incremented in real time, it presented an incremental 
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Figure 9A 
Pilot Graphics Display Sections 
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increase in each displayed state time history. The time 
histories would progress in a left-to-right movement and 
upon reaching the right most limit of the VARIABLE DYNAMICS 
section were erased from the screen and restarted at the 
left side of the electronic strip chart. Execution of a 
maneuver of the vehicle would immediately be displayed in 
the time histories displays. 
b. Sea Track/Emergency Performance Index 

The sea track portion of this display area was 
a navigational plot of the ship trajectory (Xo ¥0) (See 
Fig. 9B). This track was updated every 2 seconds as imple- 
mented in the XDS-9300 multiplex algorithm. The track 
would erase itself and begin anew at the craft's current 
position every 3.5 minutes. 

The performance of the pilot during an acci- 
dental failure situation was measured and graded through 
(TRACK DEVIATION) , 


1 
(OVERALL PERFORMANCE). All 


the use of three performance indices J 


eee cONTROL EFFORT), and J 


2 3 
performance tests were initialized such that the vehicle 


was on a heading of 000° (NORTH) and on the Xvav axis with 


a steady state velocity. 

Jas TRACK DEVIATION, was an "“integral-error- 
Squared-plus-time" performance index which penalized the 
distance that the pilot allowed the vehicle's c.g. to 


drift away from the X axis in addition to the length 
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At the systems entry into an accidental failure mode of 


operation the XDS-9300 computed Jy as 
Ceinal 5 
oe (Y_)°“ dt + TIME 
x O 
O 


Jor CONTROL EFFORT, was an index used to 
penalize the pilot for excessive control effort during the 
man controlled recovery maneuver segment of the accidental 
failure test. The Jo index was an integral "error-squared" 
index, as 


aera 
ee (oe ale 


Note that the 6 under the integral Sign is the 6 associated 
with the nozzle deflection of the three operational thrusters 
of the SES, and does not include the nozzle angle of the 
failed thruster (1.e., starboard thruster). A graphical 


interpretation of J, 1s shown in Fig. ll. 


2 
Jar OVERALL PERFORMANCE, was a magnitude 


representing the sum of Jy and J. 
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At the completion of the test the performance algorithm 

in the XDS-9300 was programmed to display the numeric 

scores representing Jy and J, and the word evaluation, 
EXCELLENT, ABOVE-AVERAGE, or GOOD below the overall performance 
heading shown in Fig. 9B. The magnitude range of the 
performance index J 3 for each of the 3 word evaluations are 


shown in Table III. 


TABLE III. OVERALL PERFORMANCE 
INDEX RATINGS 


RATING EXCELLENT ABOVE AVERAGE GOOD 


RANGE J, < 250 250) < Teele 500 < J 


3 5 3 


c. Craft Perspective/North Road/NAV Data 
Real time visual cues to the pilot were generated 

on the lower half of the pilot's graphics display console 
in the form of a perspective view of a simplified shape of 
the craft from the rear looking down and forward onto the 
craft and the area immediately ahead of the craft's path. 
An electronic north-heading road shown forward of the craft 
provided information to assist the operator in proper 
control of the vehicle. Numeric navigation data concerning 
the craft's position and the status of the nozzle angles 
and effector thrust settings were presented in the area 
immediately below the craft perspective (See Fig. 9A and 


Pag. 9B). 
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The craft perspective view presented real time 
graphical information on the vehicle's turning rate (r), 
velocity (u), pitch angle (6), and roll angle (¢). The 
electronic north road's centerline was the positive Xuay 
axis (See Fig. 1). The north road was displayed on the 
pilot's screen whenever the pilot had maneuvered the craft 
to within 1500 feet of the Xray axis and the heading of the 
vehicle was on a northerly course. At other times the north 
road was not visible to the pilot. When displayed, the 
north road provided the pilot with visual cues on the 
vehicle's heading and navigation position with respect to 
the fixed Xvay axis. When the pilot maneuvered the craft 
such that the north road was not visible graphical turn rate 
information was determined from a zero turn rate reference 
line located just above the bow of the SES. The pilot 
established a zero turn rate whenever the center point of 
the craft's bow was lined up just below the zero turn rate 
reference line. Variable speed lines moving from the 
craft's bow toward the stern gave a graphic representation 
of the forward velocity u (See Fig. 12). A sequence of 
displays during a pilot commanded control effort of left 
nozzle deflection and then right nozzle deflection is shown 
Mmeeig. 13, Fig. 14, and Fig. 15. 

The north road positioning on the pilot's 


graphic display console was solved using the following 


assumptions: 
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NOTE: 1) Craft has been commanded by pilot into a left turn 

by deflecting nozzles .32 radians (+18°). Vehicle 
is entering left turn and is passing through a heading 
of 355° with a turn rate of -1°/sec (see NAVIGATION 
DATA readout). 

2) SEA TRACK and electronic road show vehicle has dis- 
placed slightly left of north road centerline. 

3) VARIABLE DYNAMICS show velocity and acceleration time 
histories being generated as manuever develops. 


Figure 13 
SES MANUEVER (TI) 
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NOTE: 1) Pilot has held left nozzle deflection angle of .32 
rad (+18°) and vehicle has displaced left of north 
road. Vehicle has heading of 320° or 40° off the 
north road centerline. 

2) SEA TRACK confirms trajectory of craft in left turn. 

3) VARIABLE DYNAMICS show tu,tv <9, and q in steady 
State while roll rate p is still seeking steady 
state condition. 


Figure 14 
SES Manuever (IT) 
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VARIABLE DYNAMICS SEA TRACK 


U (HEAVY) | 
SeoT( LIGHT) ->—~ 

V (HEAVY) |Po)\ | 

VDOT (LIGHT) We | 

R (HEAVY) Ve | 
af OT IO OO 


SEPGHEAVY) | ,W. an, 
PDOT( LIGHT) R77 UO" 


Q (HEAVY) |. 
QDOT( LIGHT) 


SEW GEN Ci eer wnORMAN CE iN DEX 
TRACK CONTROL OVERALL 
BEY SATEC SEE PORT PERFORMANCE 


CONTROL INPUTS TIME NAVIGATION DATA 
THRUST RUDDER DELT, ELAPSED PITCH HEADING/RATE DRIFT ROLL SPEED 
400000 ae .084 Baa ro U3071 3c —s0 =o 50 
THRUST] THRUST2 THRUST3 THRUST4 RUDL RUD2 RUD 3 RUD4 
100000 ae =i =r -.95 


NOTE: 1) Vehicle has been manuevered into right turn. 
2) Vehicle heading 030° on path which will intersect 
MeteehnpOad. 


Figure 15 
SES Manuever (III) 
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1) The view of the pilot was restricted to plus and 
minus 60 degrees of the heading of the craft. This 
corresponded to the far right and far left edges 
of the pilot's display screen, respectively. 
2) The starting edge of the road would remain fixed at 
a point ahead of the vehicle whenever the vehicle was 
positioned such that the pilot could see the road in 
his restricted viewing arc. 
3) The road would give visual cues on the craft's left- 
right displacement from the road centerline and on 
the relationship between the north heading of the 
road and the heading of the vehicle. 
4) The road's shape as it appeared to the pilot would be 
a wide path terminating at an apex point in the distance. 
A non linear relationship was used to compute the 
position of the north road (See Fig. 16A and Fig. 16B). The 
Y displacement of points A, B, and C on the graphics 


NAV 

display console was made a function of Me only. The Yay 
displacement of point D on the graphics display console was 
made a function of Y 5 and the arc tangent of the craft's 
heading. The Xvay displacement of the points A, B, C, and 
D were fixed at a constant value, thus recalling Fig. 7B, 
the computation of the coordinates of the points A, B, C, 


and D was 
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Figure 16A 
North Road(I) 
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Figure 16B 
North Road(ITI) 
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Xn = mw. + X,) Yn = +.01 

x, = “Yo YR ie ts 
Xa = Bh = - X,) Yo = +.01 

X = —-(Y + 10000(tan™* Wome = 4: 
D O D : 


The displacement of points A, B, and C due to 


vehicle heading were neglected. 


E. OPERATIONAL MODES 

The RTS5D operated in any combination of modes as 
designated by the operator (See Fig. 17). Mode selection 
was accomplished via "on-off" switches located at the 
console of the XDS-9300 and at the thruster console located 
at the pilot's seat (See Fig. 5A and Fig. 5B). The RTS5D 
ran in any combination of the modes available with the 
"RESTART" mode taking priority over all other modes selected. 

The FIXED THRUST and FIXED EFFECTOR ANGLE removed the 
appropriate man-generated control and utilized a constant 
value for T and 6 respectively. 

The PRINT mode allowed variables to be printed on the 
XDS~-9300 line printer at intervals selected by the operator. 
The HOLD mode put the program into a loop whereby the 
display screens would continue to present the last image 
computed until such time as the mode was deselected. At 


this time the program would continue. 
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MODE 


10 


DESCRIPTION 


THRUST MAGNITUDE FIXED 

THRUST MAGNITUDE VARIABLE 

EFFECTOR DEFLECTION ANGLE FIXED 
EFFECTOR DEFLECTION ANGLE VARIABLE 
ACCIDENTAL FAILURE THRUST LOSS 
ACCIDENTAL FAILURE EFFECTOR REVERSAL 
PRINT 

HOLD 

RESTART ADJUST CONSTANTS IN PROGRAM 


RESTART RETAIN CONSTANTS 


Figure 17 


RTSSD Modes of Operation 
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The RESTART mode re-initialized the program and then 
provided an opportunity to the operator to change the 
magnitvde of any constant in the program. This feature 
allows the operator to carry out program modeling adjust- 
ments without the d2lays associated with re-compiling a 
modified digital program. Adjustments to constant magni- 
tudes are made via a peripheral teletypewriter (TTY-35) 
(See Fig. 5A). 

The ACCIDENTAL FAILURE FFFECTOR REVERSAL mode deflected 


the nozzle angle of the starboard thruster 4 (See Fig. 3) 


10 
to 180°. This simulated the accidental reversal of an 
peeboard thruster at full power. Simultaneously the 

program computed and displayed the performance indices 
discussed in the Software section of this report. Exit 

from this mode was accomplished when the pilot had maneuvered 
the craft back on track. At this time the program would 
"hold" while displaying the last computed values of all 
variables. The hold condition would continue until a 

restart was accomplished by the operator. 

The ACCIDENTAL FAILURE THRUST LOSS mode was entered by 
decreasing the thruster setting of any one, or combination 
of, thrusters (See Fig. 5B). Exit from this mode was 
identical to the accidental failure effector reversal 
described previously. 


The FIXED STEP SIZE INTEGRATION mode removes the real 


time clock input and utilizes an operator specified time 


no 





increment on each iteration of the program loop. MThis is 
a non real time mode. The output displays and other modes 


are functional in this mode. 
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V. RTSSD RESPONSE CHARACTERISTICS 


A. VALIDATION 

The DBSIM5D program was utilized as the benchmark to 
validate the RTS5D simulation. A DBSIM5D sequence of runs 
was made with step effector deflection angles of 5°, 10°, 
mmeeezo-, 25°, and 30°. The runs were initiated with the 
craft at 56 knts and traveling straight ahead just prior 
to effector deflection. The first peak overshoot value 
and the steady state value attained by the variables 
u (surge), v (sway), r (yaw rate), 6 (roll), and @ (pitch) 
were meed to document error in the RTS5D. The results of 
the comparison are plotted as error versus effector deflection 


Smee in Fig. 18, Fig. 19, Fig. 20, Fig. 21, and Fig. 22. 


B. ACCIDENTAL FAILURE ANALYSIS 

The response of the RTS5D when subjected to accidental 
failure mode scenarios was examined. Two primary situations 
involving the requirement for the pilot to react to the 
failure and maneuver the craft back onto itS original track 
were analyzed. The first accidental failure situation 
consisted of reversing the starboard thruster T19 (See 


frag. 3), such that 6 equalled 180° while the vehicle was 


10 
tracking straight ahead at 56 knts. The second situation 
consisted of again reversing the Tio thruster to $10 equal 


180° with the additional requirement that the failure occur 
while the vehicle was established in a steady state full right 


fete 1 (5. 7=5 5=59=65 1 )=-30°) - To assume that the "worst case"accidental 
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Figure 20 


Effector Angle vs Yaw Rate Error 
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Figure 22 
Effector Angle vs Pitch Error 
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failure scenario prevailed throughout the conduct of the 
two failure situations described above, the pilot was not 
allowed to reduce the thrust magnitudes of either the 
reversed effector Tio or the remaining thrusters To, Te. 
and To: The pilot was required to recognize the failure had 
occurred, respond with direct real time man generated control 
efforts via the "joystick" (See Fig. 5A) and maneuver the 
vehicle back onto the "north road". 

In each case the exercise was started by initializing 
the craft at the steady state conditions for either a 
Straaght ahead run or a -30° effector angle right turn. 


The appropriate steady state values for each initialization 


used are shown in Table IV. 


TABLE IV. STEADY STATE CRAFT CONDITIONS PRIOR 
TO ACCIDENTAL FAILURE TEST 


VARIABLE STRAIGHT AHEAD REGHI TURN 
u 87. ft/sec 73.99 ft/sec 
V 0. ft/sec -7.58 ft/sec 
R 0. deg/sec 2.9 deg/sec 
> 0. deg ee. deg 
0 1. deg 1.19 deg 
W QO. deg QO. deg 
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Termination of the exercise occurred whenever the pilot 
Simultaneously satisfied two recovery requirements: 
1) The craft's navigation coordinates must be within 


50 ft of the north road centerline (+X. axis). 


AV 
2) The translation velocity of the vehicle in the Y 


NAV 
direction must be less than plus or minus 1 ft/sec. 
It was considered that the craft was completely under control 
once the vehicle had fulfilled the above requirements. It 
is instructional to note that the recovery criteria not 
only required that the pilot counter any asymmetric turning 
moments and sway velocities induced by the failure, but 
that he maneuver the vehicle back to its original track. 
The results of a series of accidental thrust reversal tests 
are shown in Fig. 23A and Fig. 23B. The maximum Yxav (See 
Fig. 1) translation was 82.2 ft for the straight ahead run 
case and 571.2 ft for the right turn effector reversal case. 
It was observed that in both cases the craft retained the 
Capability to maneuver back to the original track. Further, 
the corrective action required in both cases was simply to 
command the three operational effectors to the full left 
meo-) position and hold that until the approach to the 
Original track, at which time the effectors were modulated. 
Consequently the fastest recovery can be expected to come 
from the operator who is quickest to recognize that an 


effector reversal has occurred. The initial optimal control 


algorithm in both cases was full left rudder. 
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Repeated runs made with the same operator demonstrated 
a very steep learning curve. Rapid improvement in initial 


scores was demonstrated in every case. 
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VI. CONCLUSIONS 


The ability to implement a real-time, man-controlled 
Simulator exhibiting the steady state dynamics of the 3K-SES 
has been established. By requiring the right hand side of 
the simplified equations of motion to be restricted to a 
few significant terms, that portion of the iteration loop 
delay dedicated to computation and the AD/DA process has 
been reduced to an acceptable magnitude of approximately 
35 msecs. The remaining portion of the total iteration time 
of approximately 84 msec was required for the subroutine's 
outputting the "alpha-numeric" presentation of the craft's 
States onto a graphics display console. This low value of 
graphics presentation time resulted from the utilization 
of time multiplexing of the data to be displayed. Using 
100 msecs aS a maximum allowable iteration time additional 
terms for the right hand side can be incorporated to generate 
amore accurate simulation. From Section III on rectuangular 
integration error analysis it was concluded that the RTS5D 
loop delay must remain less than 100 msec to remain credible. 
This concurs with those conclusions reached by A.L. Greensite 
in Ref. 9. Once the 100 msec loop delay point has been 
reached further expansion of the right hand side of the equa- 
tions of motion would come from development of a more 
efficient graphics system. It was observed that while the 


graphics console utilization required a high cost in terms 
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of loop delay, the interface of realistic man-generated 
control efforts could be accomplished in no other available 
manner due to limitations on existing facilities. 

Results of operational use of the RTS5D indicate the 
craft perspective and north road algorithm developed in 
Section IV-D and used to generate real time visual cues 
for the pilot to be extremely effective in displaying craft 
heading changes and vehicle displacements in Yo: 

The results in Section V-A have shown the RTS5D to be 
within acceptable error limitations of less than 2% on the 
average when compared to the DBSIMS5D program in the moderate 
speed range of approximately 56 knts and over the full 
effector deflection range of plus or minus 30 degrees. 

The results in Section V-B have shown that the 3K-SES 
is capable of being controlled satisfactorily by man during 
accidental failures involving both loss of thrust and 
effector reversals. From these results it was shown that 
the optimal control effort during the initial stages of an 
effector reversal failure was simply to apply full opposing 
effector angle on the remaining operational effectors. This 
input was held by the pilot until the craft returned to the 
vicinity of the original track at which time the operational 
effectors deflection angles had to be man-manipulated 
Slightly to maintain track. Consequently computer assisted 


maneuvering controls would be somewhat helpful but do not 
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appear to be a hard requirement during effector reversal 
scenarios in the 3K-SES. 

It was also concluded that the RTS5D program can be 
readily expanded to a full 6 degree of freedom simulator 
with the addition of simplified equations representing the 


heave motion characteristics of the 3K-SES. 
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VII. RECOMMENDATIONS 


The following recommendations are made for future 
Studies utilizing the RTS5D to simulate the 3K-SES: 

Weerne Current right hand side of the craft's equations 
of motion require expansion to include more terms to 
accurately represent the 3K-SES over a wider operating range. 

2. Incorporate the sixth degree of freedom, heave, 
into the simulation. 

3. Include "sea-state" input whenever more computational 
time becomes available. 

4. Program the propulsion plant dynamics to allow 
constant power aS an input in lieu of constant thrust. 

5. Design a dedicated micro-processor or mini-computer 
to assist the operator in maneuvering control with the intent 
of more fully examining the computer assisted controls 
requirement in the 3K-SES. 

6. Incorporate actual 3K-SES control hardware into the 
RTS5D to allow more realistic "hands on" testing of not only 
the control hardware itself but also the operator personnel. 
Currently a project is being designed to test operator 
response under various command and failure conditions. 


Results will be reported in the near future. 
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All 
AZ1 
A22 
A31 
A33 
A34 
A4l 
A51 
A6l 
A7l 
A81 
A91 
AASTOP 
ABSYD 
ABSYOD 
AD 
A1COA 
AM 

f LTCH 
APOTK5 
APRINT 


AR 


APPENDIX A 


RTSSD PROGRAM NOMENCLATURE 


A/D trunk 500 voltage 

A/D trunk 501 voltage 

Added mass coefficient in yaw moment equation 
A/D trunk 502 voltage 

Added mass coefficient in roll moment equation 
Added mass coefficient in pitch moment equation 
A/D trunk 503 voltage 

A/D trunk 504 voltage 

A/D trunk 505 voltage 

A/D trunk 506 voltage 

A/D trunk 507 voltage 

Pye trunk 51O voltage 

performance index limit 

IY, | 

Fa 

array of A/D lines 

speed line index counter 

mass 

puecheangle in degrees 

weltageron pot E of thruster console 


# iterations between print execute commands 


analog turn rate voltage 
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ARATE 
ARDOT 
AROLL 
ASTOP 
AU 

AUDOT 
AV 


AVDOT 


AXST 
AY 
AYD 
AYST 
B101 
Bid 1 
ed 1 
pe 1 
ES 1 
B41 
B51 
ED 
Bol 
B/1 
B81 
eo 1 


BETA 


turn rate in degrees 
analog r voltage 


roll angle 


performance index J 


1 ‘sedate 


analog velocity voltage 


analog u voltage 


analog sway velocity voltage 


analog v voltage 


X-Y plotter scale factor 


KY plotter X, coordinate 


initial condition X-coordinate on sea track 


mY OlOtter Scale factor 


MY Totter Yo coordinate 


initial condition Y-coordinate on sea track 


D/A 
D/A 
D/A 
D/A 
D/A 
D/A 
D/A 
D/A 
D/A 
D/A 
D/A 


D/A 


trunk 
Eaunk 
cesta a 
ten kK 
trunk 
trunk 
trunk 
trunk 
trunk 
i eunrnue 
eran kK 


trunk 


T-432 
T-420 
T-433 
T-421 
eee 
T-423 
T-424 
Tea 
T-426 
T-427 
T-430 


Raa 31 


drift angle 
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BOAT 
CDP 

CDX 

CDY 
CDZP 
DAL 
DELT 
DELTA 
DTIME 
DTIMPLT 
DXPLOT 
DYPLOT 
DYREPET 
EFLG 
ETIME 
EYREPET 
EXPLOT 
EYPLOT 
Fl 

F2 

F3 

FSP 

FSS 
FTIME 
FTIMPLT 
FXPLOT 


pe PLOT 


Subroutine to generate perspective of SES 
bow real forces lumped coefficient 
drag forces lumped coefficient 
sway forces lumped coefficient 
Sidewall rolling moment lumped coefficient 
digital to analog call 

iteration loop time 

commanded iteration loop time 

v time 

scaled value for v time history 
same as DTIMPLT 

scaled value for v 

restart value for v 

Sruonw ctlag for TEXTO call 

r time 

restart value for r 

scaled r time 

scaled r 

stern buoyancy force 

bow seal pitch force 

bow buoyancy force 

port sidewall buoyancy force 
starboard sidewall buoyancy force 
7 time 

scaled v time 

same as FTIMPLT 


scaled v 
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FYREPET 
G 
GRAPHO 
GTIME 
GTIMPLT 
GXPLOT 
GYPLOT 
GYREPET 
HEAD 
HOLD 
HTIME 
HTIMPLT 
HXPLOT 
HYBETA 
HYBRID 
HY PLOT 
mer ST 
HYR 

HY RDOT 
HYREPET 
HYSESX 
HYSESY 
HYU 
HYUDOT 
HYV 


HYVDOT 


restart value for v 


gravitational acceleration 


subroutine to project image 


R time 


scaled 


(non alpha-numeric) 


R time 


same as GTIMPLT 


scaled 


R 


restart R value 


craft heading 


subroutine to freeze display/program 


u time 


scaled 


u time 


same as HTIMPLT 


analog 
analog 
analog 
analog 
analog 
analog 
analog 
analog 
analog 
analog 
analog 
analog 


analog 


B 

Plraqme@eetrvcient (0 or i) 

y axis sea track coordinate 
(radians) 


heading 


turn rate (rad) 


a 
restart y coordinate on sea track 
Ver ODOsLeLON Of CHabt 


MepOStELON. Of Claim 


u 
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HYXO analog Xo 

HYXODOT analog Mie 

HYYO analog YG 

yy ODOT analog Y, 

EGO multiplex graphics index 
ICOA speed line index 

ICRAFT vehicle designator 

IER error flag 

IIPR time history limit 

IIPRD time history limit 

IIPRE time history limit 

ITIPRF time history limit 

IIPRG tame History limit 

IIPRH time history limit 

IIPRO time history limit 

IIPRP time history limit 

IIPRQ time history limit 

IIPRR time history limit 

IIPRS time history limit 

EX Prat ecounter 

IJ craft perspective counter 
IJA craft perspective counter 
IJAD time history counter 
IJAE time history counter 
IJAF time history counter 
IJAG time history counter 
IJAH time history counter 
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IJAO 
IJAP 
IJAQ 
IJAR 
IJAS 
IJD 
IJE 
GIF 
IJG 
IJH 
IJO 
LIP 
IJQ 
IJR 
IJS 
IL 
ILA] 
ILA2 
ILA3 
ILA4 
ILA5 
ILA6 
ILA7 
ILA8 
ILA9 


ILA1LO 


time 
time 
time 
time 
time 
time 
time 
time 
time 
time 
time 
time 
time 
time 


time 


history 
history 
ls EO ry, 
history 
history 
history 
history 
history 
history 
history 
history 
history 
history 
history 


history 


Counter 


counter 


counter 


counter 


counter 


counter 


counter 


counter 


counter 


counter 


counter 


counter 


counter 


counter 


counter 


graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 


graphics 


console 
console 
console 
console 
console 
console 
console 
console 
console 
console 


console 


display 
display 
display 
display 
display 
display 
display 
display 
display 
display 


display 


G2 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 





ILAI1 
ILA12 
ILA13 
ILA14 
ILAL5 
ILA16 
ILA1?7 
ILA18 
ILA19 
ILA20 
ILA21 
ILA22 
ILA23 
ILA24 
ILA25 
ILA26 
ILA27 
ILA28 
ILA29 
ILAP 
IN 


IPACK 


ei RET 


Pe LOT 


eR 


ee RD 


graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 


graphics 


console 
console 
console 
console 
console 
console 
console 
console 
console 
console 
console 
console 
console 
console 
console 
console 
console 
console 
console 


console 


display 
display 
display 
display 
display 
display 
display 
display 
display 
display 
display 
display 
display 
display 
Agaclasy 
display 
display 
display 
display 


display 


failure mode counter 


subroutine to load display arrays 


performance index 


array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
SaReny 


array 


graphics console display array 


time history counter limit 


time history counter limit 
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[PRE 
ie RE 
IPRG 
IP RH 
[PRO 
ire RP 
me RO 
IPRR 
fe RS 
ISWA 
oe XT 
IVEIWA 
IVIEWAA 
IVIEWB 
IVIEWC 
IVIEWD 
IVIEWE 
IV TEWF 
IVIEWG 
IV IEWH 
IVIEWO 
IVIEWP 
IVIEWOQ 
IVIEWR 
IVIEWS 


IVIEWZ 


time 
time 
time 
time 
time 
time 
time 
time 


time 


display array 
Subroutine to display alpha-numeric data 


craft perspective/north road array 


NiUsSteLy. 
history 
history 
pau cte ny) 
history 
history 
history 
history 


history 


counter 


COUNnter 


COuUnEcerT 


counter 


COuUnRECY 


counter 


counter 


counter 


counter 


sea track array 


sea track axis 


array 


u time 
v time 
mei me 


time 


Ge 


time 


Ke 


time 


ce 


p time 
g time 
p time 
q time 


analog 


history 
history 
history 
history 
history 
history 
ES EOny 
history 
his cory 


history 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


mast 
aie 
Por eete 
ae m ante 
limit 
eee 
Tana 
Meera 


deamane 


sea track array 


time history axis array 
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IX 


JD 
JE 
om 
JG 
JH 
JO 
oP 
JQ 
JR 
JS 
L3 
LD 
or 


LZ 


NAD 

NDA 

OO 
OTIME 
OTIMPLT 
exPLOT 
OYPLOT 
wrnePET 
E 


PB 


X AaxX1S moment of inertia 


time 
time 
time 
time 
time 
time 
time 
time 
time 
time 


time 


history 
Risto ry 
history 
history 
hus tony 
history 
history 
history 
history 
Huston y 


history 


counter 


COumee th 


counter 


counter 


counter 


COuntCer 


counter 


counter 


counter 


counter 


counter 


lever arm for buoyancy force 
draft of sidewall 

Side thrust component lever arm 
lever arm for sway drag force 
Slockt interrupt count 

AD trunk line array 

DA trunk line array 

lever arm for thrust side component 
p time 

scaled p time 

same as OTIMPLT 

scaled Dp 

restart p value 

eM IE verelaaxe! 


plenum pressure 
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PDOT 
PDOTM 
PHI 
PMAX 
POTK 
POTS 
POTSA 
rogk 

eri ME 
Eee P TT 
EAP LOT 
ere LOT 
Pere ePET 
Q 

QDOT 
QDOTM 
QMAX 
QTIME 
QTIMPLT 
OxP LOT 
oye LOT 
OL REPET 
R 

RDOT 


RDOTM 


roll acceleration 
max roll acceleration 
roll angle 

max roll rate 
thruster console array 
Control stick array 
function box array 
craft heading 

g time 

scaled q time 

same as PTIMPLT 
scaled q 

restart gq value 
patch rate 

pitch acceleration 
max q 

max q 

p time 

scaled p time 

same as QTIMPLT 
scaled p value 
restart p value 
Geaet Lurbn rate © 


R 


max 


READCLOCK subroutine to sample real time clock 
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RESET 
RHO 
RMAX 


ReERE 1 


RPERF2 
PEER 3 
RTIME 
Pao P LT 
RX 

Rae LOT 
RY 

Ry PLOT 
RYREPET 
RZ 

Sl 

S2 

S3 

S4 
SESX 
SESY 
lal 
SF14 
SP 
SPDLIN 
SPEED 


STARTCLOCK 


subroutine to reset analog computer 
density of water 


maximum yaw rate of craft 


q time 
scaled q time 


Ty 


same as RTIMPLT 


ay 


scaled gq value 
restart gq value 


Ty 


Ty lever arm 


Ti 9 lever arm 


Ts lever arm 


To lever arm 
scaled X 

O 
scaled Y 

O 
scale factor in time history arrays 
scale factor in time history arrays 
Beta influence coefficient 
speedline lower screen limit 
craft's velocity in knts 


Subroutine to start clock 
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STIME 
saMPLT 
eee LOT 
exPLOT 
wene. PET 
Sen EPET 
a 

£10 
ELOC 
TLOMAX 
LOS 

7 

7c 

T 7MAX 
iS 

28 

T8C 
T8MAX 
T8S 

9 

IC 
TIMAX 
oS 
TEXTO 


TFORW 


analog time 

scaled analog time 

HYSESX 

neobo Y 

restart value for HYSFSX 
Bestarnt value for HYSESY 

total thrust 

thruster number 4 

forward vector component of T1LO 
max T1lO value 

Side vector component of T10 
thruster number 1 

forward vector component of T7 
max T7 value 

Side vector component of T7 
thruster number 2 

forward vector components of T8 
max value of T8 

Side vector component of T8 
thruster number 3 

forward vector component of T9 
max value of T9 

side vector component of T9 
subroutine to display alpha-numeric data 


summation of forward thrust components 
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PAHETA 
iM 
TIME 
fee? LT 
TINT 
feet LE 
Pet LEO 
PITLEL 
TMAX 
TSIDE 
TYAW 

U 

UDOT 
BDOTM 
UMAX 


UMUL 


VCD 
VDOT 
VDOTM 
VMAX 


VS 


WDOT 
WE 


WRITECLOCK 


Prech angie 

u time 

time 

scaled TIM 

real clock interrupt frequency 
program name display array 
program name display array 
program name display array 
max T 

sum of side thrust components 
sum of yaw moment 

forward velocity (surge) 

u 

max u 

max u 

scale factor in speedline generation 
side velocity (sway) 

thruster console pot array 

Vv 

max v 

max v 

total velocity 

heave velocity 

heave acceleration 

width of bow seal 


subroutine to assign clock interrupts value to 
variable 


oe 





X 

XO 
XODOT 
XDX 
XGX 
XPLOT 
XREPET 
XST 
YO 
mODOT 
YHIGH 
we LOT 
ie ET 
YST 

Z 

Z10 
Z2 

Z3 

Z4 

Z5 

Z7 
Z/71I 
Z8 
Z8I 
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lever arm for sway drag forces 


boat perspective scale factor 
boat perspective scale factor 
display array 

reset value for SESX 

scale factor for analog plot 


Yo 
EG 

point D vertical position constant 
display array 

reset value for SESY 

scale factor for analog plot 
operational effector angle 
effector number four angle 

upper effector angle limitation 
lower effector angle limitation 
zero angle upper limit 

zero angle lower limit 

effector number one angle 

modified effector number one angle 
effector number two angle 


modified effector number two angle 


effector number three angle 
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Z9I modified effector number three angle 


ZAB scale factor for 210 
ZI modified Z 
ZZZZ fixed step effector angle input 
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